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(Dated: April 9, 2015) 

Unpolarized and beam-polarized four-fold cross sections e'p'^ reaction 

were measured using the CLAS detector and the 5.75-GeV polarized electron beam of the Jefferson 
Lab accelerator, for 110 {Q ^bins over the widest phase space ever explored in the valence- 
quark region. Several models of Generalized Parton Distributions (GPDs) describe the data well at 
most of our kinematics. This increases our confidence that we understand the GPD H, expected to 
be the dominant contributor to these observables. Through a leading-twist extraction of Compton 
Form Factors, these results reveal a tomographic image of the nucleon. 

PACS numbers: 12.38.-t, 13.40.Gp, 13.60.Fz, 13.60.Hb, 14.20.Dh, 24.85.-|-p 


The internal structure and dynamics of the proton, the 
nucleus of the most abundant chemical element in the 
visible universe, still remain a mystery in many respects, 
more than 40 years after the evidence for its quark and 
gluon substructure. How are the spatial and momentum 
distributions of the quarks and gluons (Le., the partons) 
correlated inside the nucleon? How do the partons con¬ 
tribute to the bulk properties of the proton (mass, spin, 
charge,...)? These are some fundamental questions at the 
intersection of nuclear and particle physics which are still 
to be resolved. 

In order to tackle these essential issues, a large ex¬ 
perimental program was launched worldwide at Jeffer¬ 
son Lab (JLab), COMPASS and HERA, facilities using 
multi-GeV electromagnetic probes, to study deeply vir¬ 
tual Compton scattering (DVCS). In the valence-quark 
region, this corresponds to Compton scattering at the 
quark level, with the incoming photon radiated from 
the lepton beam. As in the study of atomic or nu¬ 
clear structure, the energy and angular distributions of 
the scattered photon reflect the distribution in momen¬ 
tum and/or space of the targets, which in our case are 
the quarks inside the proton. At JLab, electron beams 
are used and the reaction to study proton strncture is 
ep —>■ e'p'y. It was shown w that this process, at 
sufficiently large squared electron momentum transfer 
= — (e — e')^ and small squared proton momentum 
transfer t = {p — p'Y terms of the electron and pro¬ 
ton four-vectors), could be interpreted in the framework 
of Quantum Chromodynamics (QCD), the fundamental 
theory governing the interaction of quarks and gluons, as 
the product of the elementary Compton scattering at the 
quark level 'y*q —>■ 7(7 with factorizable structure func¬ 
tions called Generalized Parton Distributions (GPDs). 

Figure (left) illustrates GPD QGD factorization for 
the DVCS process. In a frame where the nucleon moves 
at the speed of light in a given direction, a quark with 
longitudinal momentum fraction a; -I- ^ absorbs the vir¬ 
tual photon and, after radiating the final-state photon. 
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cess at sufficiently large and small |t|, as predicted by 
the QCD factorization theorem. Right: the Bethe-Heitler 
(BH) process, which interferes with the DVCS process in the 
ep —^ e'p^7 reaction. 

the same quark returns into the nucleon with a longi¬ 
tudinal momentum fraction x — plus some transverse 
kick included in t. The GPDs are functions of x, and 
t, and represent the probability amplitude of such a pro¬ 
cess. The variable ^ is related to the Bjorken variable xb- 
^ ~ 2-xb ’ Xb — 2 ^ with the proton mass M and 

V = Ef. — Ef.'. Thus, it is determined by the scattered- 
electron kinematics. The quantity x is not measurable 
in the DVGS process. At leading-order QCD, GPDs do 
not depend on . At leading-twist QCD, four GPDs 
enter the description of the DVCS process: H, H, E and 
A, representing the four independent helicity-spin transi¬ 
tions of the quark-nucleon system between the initial and 
final states. The GPDs are QGD matrix elements that 
integrate and project on a few variables the full com¬ 
plexity of the quarks’ and gluons’ dynamics within the 
nucleon. 

The GPDs embody the longitudinal momentum distri¬ 
bution of the quarks in the nucleon, their transverse spa¬ 
tial distribution, and the correlation between these two 
distributions. One uses the term nucleon tomography as 
one can probe the transverse size of the nucleon for dif¬ 
ferent quark longitudinal-momentum slices. For details 
on the GPD formalism, see the reviews [6HII]. 

In the ep —>■ eVb reaction, the DVGS process inter¬ 
feres with the well-known Bethe-Heitler (BH) process, 
illustrated in Fig. (right), where the final-state photon 
is radiated by the incoming or scattered electron. 

Extracting the GPDs from the DVCS process requires 









3 


measuring a series of observables for the ep —>■ e'p'j reac¬ 
tion over the broadest kinematic domain possible. Sev¬ 
eral observables, such as the unpolarized cross section 
and polarized beam or/and target asymmetries, are nec¬ 
essary in order to separate the four GPDs. Each observ¬ 
able is sensitive to a particular combination of GPDs. 

This paper presents a major contribution to this global 
and long-term endeavour: the extraction of the ep —>■ 
eV7 DVGS-I-BH) unpolarized and beam-polarized 
cross sections over the widest phase space ever explored 
in the valence-quark region, with 110 bins 

covering: 1.0 < < 4.6 GeV^, 0.10 < xb < 0.58, 

and 0.09 < —t < 0.52 GeV^. In this kinematic do¬ 
main, our results strongly enhance the existing set of 
measurements of the ep —> e'p'^ reaction which con¬ 
sists of: four {Q‘^,XB,t) bins of unpolarized cross sec¬ 
tions and 12 bins of beam-polarized cross sections mea¬ 
sured by the JLab Hall A collaboration [12], 57 bins of 
beam-spin asymmetries measured by the CLAS collabo¬ 
ration m, and 18 bins of longitudinal target- and beam- 
target double-spin asymmetries measured by the GLAS 
collaboration dSds] (in addition to the handful of GLAS 
pioneering data points of Refs [TMT 8 ] i. 

The experiment took place at JLab during three 
months in 2005, using the 5.75-GeV polarized elec¬ 
tron beam (79.4% polarization), a 2.5-cm-long liquid- 
hydrogen target, and the Hall B large-acceptance 
CLAS spectrometer m, operating at a luminosity of 
2x10^^ cm“^s“^. A specially designed electromagnetic 
calorimeter (“inner calorimeter”, IC [T3]) was added to 
the CLAS detector and allowed the detection of photons 
for polar angles from about 5° to 16°, with full azimuthal 
coverage. A solenoid magnet was installed around the 
target in order to magnetically shield the IC from the 
copious Mpller background stemming from the target. 

The first step of the data analysis was to select events 
with at least one electron, one proton, and one photon 
in the final state. Electrons were identified by signals in 
the CLAS drift chambers, scintillators, Cherenkov coun¬ 
ters, and the standard CLAS electromagnetic calorime¬ 
ters. Protons were identified by the correlation between 
their measured momentum and velocity. The highest- 
energy particle detected in the IC was considered as a 
photon candidate. Once these three final-state parti¬ 
cles were selected and their 3-momenta determined, the 
exclusivity of the ep —>■ e'p'j reaction was ensured by 
applying 3a cuts on the following four variables: the 
squared missing mass MM^,p, of the (e'p'X) system, the 
coplanarity angle A(j), i.e., the angle between the {'y*,p') 
and ( 7 *, 7 ) planes, the missing transverse momentum of 
the {e'p'^) system, and the angle O-yx between the mea¬ 
sured photon and that predicted by the kinematics of 
the {e'p'X) system. We also selected the particular kine¬ 
matics: IT > 2 GeV, where IT^ = s = ( 7 * +p)'^, to 
minimize contributions from radiative decay of baryonic 
resonances, and > 1 GeV^ in order to be in the deep 



FIG. 2. (Color online) Two of the four variables on which 
So- exclusivity cuts (vertical lines) were applied to select the 
ep —>■ e'p'y reaction: A(j) (left) and 9~fX (right). Black solid 
distributions correspond to the events with at least one elec¬ 
tron and one proton detected in CLAS, and one photon 
detected in the IC, after applying the exclusivity cuts on 
MMg/p/. Each blue shaded distribution corresponds to the 
events remaining after applying the exclusivity cuts on all 
the variables except for the plotted one. 


virtual regime. As an example, Figure shows the effect 
of two of the four exclusivity cuts. 

Under these conditions, we ended up with about 
300,000 events. Figure shows the resulting {Q'^,xb) 
and {—t,XB) kinematic coverages of the data and the 
adopted binning [21 {Q‘^,xb) bins and 6 t bins], which is 
finer than the one used in Ref. m- Note that the bins 
and results presented here are limited to the jtj region 
below 0.52 GeV^ while the actual coverage of the data 
goes beyond 1 GeV^. The ep —>■ e'p'^ cross sections vary 
very rapidly with kinematics, primarily due to the BH 
process. In order to minimize the uncertainties related 
to the knowledge of the kinematics, we chose to minimize 
the size of our bins, while keeping comparable statistics 
in each bin. 



FIG. 3. (Color online) The {Q^,xb) (left) and {—t,XB) (right) 
kinematic coverages of this experiment, with the correspond¬ 
ing binning. 


Due to the azimuthal symmetry when using an unpo¬ 
larized target, the ep e'p '7 reaction depends on four 
independent variables. For the study of GPDs, the most 
appropriate ones are xb, t and 4>, where (j) is the 
azimuthal angle between the (e, e') and {'y*,p') planes 
around the virtual photon direction. We have thus ex- 
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tracted four-fold cross sections as follows: 




^ep—^e'p'^ 


N, 


ep—ye'p''y 


dQ^dxsdtdcj) Cmt^Q^ b ^tAcj) AcC Frad 


( 1 ) 


In Eq. ^ Nep^e'p'-y is the number of ep e'p'q events 
in the 0) bin. The aforementioned exclusiv¬ 

ity cuts do not fully select a pure sample of DVCS-I-BH 
events. We evaluated the contamination from the ep —>■ 
e'p'ir^ channel where one photon of the 7r° decay can es¬ 
cape detection, using a combination of ep —)■ e'p'ir^ mea¬ 
surements and Monte-Carlo simulations. On average, 
this contamination is less than 9% and was subtracted 
on a bin-by-bin basis. The four-dimensional accep¬ 
tance/efficiency of the CLAS detector, Acc, for the ep —>■ 
e'p'^ reaction was determined for each (Q^, ccb, t, bin 
by generating more than 200 million DVCS-I-BH events, 
using a realistic Monte-Carlo generator. The events were 
processed through the GEANT simulation of the CLAS 
detector, and the same reconstruction and analysis codes 
that were used for the data. The event generator includes 
radiative effects so that Acc also corrects for a part of 
the real internal radiative effects. The factor Frad cor¬ 
rects, for each (Q^, ccb, t, bin, for the virtual internal 
radiative effects and the remainder of the real internal 
radiative effects, which can be both calculated theoret¬ 
ically HO]. The product (AQ^AxsAfAc/i) corresponds 
to the effective hypervolume of each bin. Finally, Cint 
is the effective integrated luminosity, corrected for the 
data acquisition dead time, which was deduced from the 
integrated charge of the beam measured by a Faraday 
cup. In addition, we applied a global renormalization 
factor of 12.3%, determined from the analysis of the elas¬ 
tic scattering ep —)■ e'p', by comparing the experimental 
cross section to the well-known theoretical one. This fac¬ 
tor compensates for various kinematic-independent inef¬ 
ficiencies, not well reproduced by the simulations. 

Figure 1^ shows, for two selected {Q‘^,xb) bins in dif¬ 
ferent parts of the phase space, the (/>-dependence of the 
ep —)■ e'p'^ unpolarized cross section and beam-polarized 
cross-section difference. The latter of these two observ¬ 
ables is defined as follows: 


A(dV) = - 




ep—^e'p''^ 




ep—^e'p'") 


dQ'^dxBdtdcj) dQ‘^dxBdtd(j) 


( 2 ) 


where the arrows correspond to beam helicity states -I- 
and —. For each of these (Q'^.xb) bins, three selected t 
bins are shown. In Fig. the black error bars show the 
statistical uncertainties of the data [13.9% on the unpo¬ 
larized cross section on average, over the 110 (Q^,XB,t) 
bins] and the blue bands show the systematic uncertain¬ 
ties [14% on the unpolarized cross section on average]. 
The contributions to the latter include the uncertain¬ 
ties on the beam energy and therefore the kinematics 
and associated corrections (5.7% on average), the accep¬ 
tance correction (5.3%), the global renormalization factor 


-t=0.153GeV2 -t=0.262GeV2 -t=0.447 GeV^ 



-t=0.204 GeV^ -t=0.262 GeV^ -t=0.448 GeV^ 



FIG. 4. (Color online) Top six plots: unpolarized cross sec- 

d‘^<7 

TOw) and beam-polarized cross-section 
difference A{d'^a) for the ep —>■ e'p'y reaction, as a function 
of <j>y for (Q^,a;s)=(1.63 GeV^, 0.185) and for 3 —t values: 
0.153, 0.262 and 0.447 GeV^. Bottom six plots: same observ¬ 
ables for (Q^ 2 ;s) = (2.78 GeV^ 0.335) and -t=0.204, 0.262 
and 0.448 GeV^. The green long-dashed curves show the BH 
contribution only. The other curves correspond to the pre¬ 
dictions of four GPD models: VGG [6l 1211 I22j (blue solid 
curves), KMS (cyan dash-dotted curves), and two ver¬ 
sions of the KM model [24l|2^, KMIO (red dotted curves) 
and KMlOa (red short-dashed curves). The blue bands show 
the systematic uncertainties. 


(5%), the exclusivity cuts (3.5%), the radiative correc¬ 
tions (2.2%), the particle selection (1.6%), and the 7r° 
background subtraction (1%). 

The unpolarized cross sections peak towards (j)=Q° and 
360° due to the BH process for which the final-state pho¬ 
ton is predominantly emitted in the direction of the initial 
or scattered electron. This is quantitatively confirmed by 
the calculations shown in Fig. where the green curves 
show the BH contribution only. The difference between 
the BH curves and the data can thus be attributed to the 
DVCS process, and therefore linked to GPDs. We dis¬ 
play in Fig. [^calculations of four GPD models, listed in 
the caption. The modeling of the GPDs in the VGG and 
KMS models is based on the Double-Distribution repre¬ 
sentation [U m HTj. The VGG calculations in Fig. |4| 
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only include the contribution of the GPD H as the in¬ 
clusion of the other GPDs barely changes the results. 
The KM model is based on the Mellin-Barnes represen¬ 
tation [231HH]- The KMIO version of the model includes 
contributions from all four GPDs for which the free pa¬ 
rameters were fitted to the JLab mils], HERMES HH] 
and ZEUS/Hl [301 EH data. In that work, it was found 
that it is possible to fit the JLab Hall A unpolarized cross 
sections only at the price of the introduction of a very 
strong H contribution [32]. The KMlOa version is based 
on a fit which excludes the JLab Hall A unpolarized cross 
sections m and sets H to zero. Note that none of these 
four models has been tuned to our data. 

Eigure [^ shows that the predictions of standard GPD 
models like VGG, KMS, and KMlOa, whose compati¬ 
bility is remarkable in spite of their different approaches, 
are in good agreement with our unpolarized cross-section 
data. In contrast, we see that the KMIO version, which 
includes the strong H contribution, tends to overestimate 
our data. Over our 110 {Q^,XB,t) bins, the average 
value per degree of freedom [33| is 1.91 for VGG, 1.85 
for KMS, 1.46 for KMlOa, and 3.94 for KMIO. We can 
therefore conclude that standard GPD models with a 
dominant contribution of the GPD H to the unpolarized 
cross section, i.e., without the introduction of a strong H 
contribution, describe the data well. Moreover, the dis¬ 
agreement between our data and the KMIO model, which 
instead matches the Hall A results, might reveal an in¬ 
consistency between the two sets of data. As a check, 
we performed a dedicated data analysis using the exact 
same {Q^,XB,t) bin limits as those used for the Hall A 
analysis (Q^=2.3 GeV^, 2:5=0.36, and —t =0.17, 0.23, 
0.28 and 0.33 GeV^). However, in this limited and par¬ 
ticular (Q‘^,XB,t) region, the comparison is hampered by 
our large statistical uncertainties and lack of ^-coverage 
around (j) = 180°. Thus no conclusion can be drawn from 
this comparison. The Hall A experiment was run at a lu¬ 
minosity almost three orders of magnitude larger than 
ours, but in a much more limited phase space. 

In general, the four models, including KMIO, give a 
good description of the beam-polarized cross-section dif¬ 
ference and the data barely allow one to distinguish one 
model from another. Over our 110 {Q^,XB,t) bins, the 
average value per degree of freedom [53] is 1.40 for 
VGG, 1.84 for KMS, 1.06 for KMlOa, and 1.20 for KMIO. 

Einally, we attempted to extract directly some GPD 
information from these two sets of observables. We used 
the local-fitting procedure developed in Refs. [33H37] . At 
leading-twist and leading-order, this procedure uses well- 
established DVCS amplitudes and does not depend on 
model parametrizations of the GPDs. We fit simultane¬ 
ously the (/)-distributions of our unpolarized and beam- 
polarized cross sections at a given {Q'^,XB,t) kinematic 


Q^=1.11GeV^ Q“=1.63GeV“ Q^=2.23 GeV“ 



0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 


-t (GeV^) -t (GeV^) -t (GeV^) 

FIG. 5. (Color online) Results of the CFF fit of our data 
for Him (upper panel) and Hue (lower panel), with only the 
GPDs H and H, for three of our (Q^,xb) bins, as a function 
of t. The blue solid curves are the predictions of the VGG 
model. The black dashed curves show the fit of the results by 
the function Ae^’*. 

point by the eight (real) quantities: 


In Eq.[^ F = H,H, E, E, the top and bottom signs apply 
to the unpolarized [H, E) and polarized (77, E) GPDs 
respectively, and V is the principal value integral. These 
quantities are called Compton Form Factors (CFFs) [33] 
in Refs. [341137] and “sub-CFFs” in Ref. [39]. The only 
model-dependent input in the procedure is that the CFFs 
are allowed to vary in a very conservative limited range, 
±5 times the CFFs from the VGG model [22]. In spite of 
the underconstrained nature of the problem, i.e., fitting 
two observables with eight free parameters, the algorithm 
manages in general to find well-defined minimizing values 
for Him and Hue- The reason is that the two observables 
that we fit are dominated by the contribution of the GPD 
77. 

Ideally, one would like to fit all CFFs. However, with 
only two observables in this case, this leads to too large 
uncertainties. We therefore present in Figure]^ for a se¬ 
lection of three of our 21 ,xb) bins, the 7-distribution 
of the fitted Him and Hue, computed nedecting the con¬ 
tributions associated with E and E. Fig.[^also shows the 
predictions of the VGG model, which overestimates the 
fitted Him at the smallest values oi xb- 

We have fitted, in Fig. [^ the 7-dependence of Him 
by the function Ae^* with the normalization A and the 
slope b as free parameters. Keeping in mind that the 
values are different for the three xb bins, the results of 
these fits show that A and b increase, in a systematic way, 
with decreasing xb- Under the hypothesis of neglecting 
higher-twist and evolution effects as well as deskew¬ 
ing effects [30] , these behaviors might reveal tomographic 
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features of the quark content of the nucleon. Under the 
mentioned conditions, b is related to the transverse size 
of the nucleon. Our data therefore suggest, over the xb 
range explored in this experiment, that the size of the 
nucleon increases as lower momentum fractions (propor¬ 
tional to Xb) are probed. The rising of A reflects the 
increase of the partonic content of the nucleon as lower 
Xb values are probed. Njie does not lend itself easily to a 
simple interpretation as it involves a weighted-integration 
of the GPD iJ over x. Nevertheless, its extraction is of 
great use to constrain models. 

In conclusion, we have measured the unpolarized and 
beam-polarized four-fold cross sections for the 

ep —>■ e’p'^f reaction, for 110 {Q^,XB,t) bins, over the 
widest phase space ever covered in the valence-quark re¬ 
gion. This large set of new data will provide stringent 
constraints on GPD models. We have shown that three 
well-known GPD models, VGG, KMS, and the KMlOa 
version of the KM model, describe the data well with¬ 
out additional inputs. The model interpretation of the 
present results favors a smaller deviation from the pure 
BH process around (j) = 180° than suggested by the Hall 
A data. Within such models, this reinforces the expec¬ 
tation of the i7-dominance in the unpolarized cross sec¬ 
tion. We have also extracted the and Hue CFFs 
from our data, fitting simultaneously both cross-section 
observables. Under some assumptions, our results tend 
to show that the nucleon size increases at lower parton- 
momentum values, thus revealing a first tomographic im¬ 
age of the nucleon. 

The full data set is available at m- A long article 
is in preparation and will include the results for all our 
{Q^,XB,t) bins and a more detailed description of the 
data analysis. 
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